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ABSTRACT 

We introduce a method for the study of CP-violating asymmetries in 
tagged states of neutral B mesons with arbitrary coherence properties. 
A set of time-dependent measurements is identified which completely 
specifies the density matrix of the initial state in a two-component 
space with basis vectors and B , and permits a determination of 
phases in the Cabibbo-Kobayashi-Maskawa matrix. For a given tag- 
ging configuration, the measurement of decays both to flavor eigen- 
states and to CP eigenstates provides the necessary information. 
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The search for CP violation in systems involving B mesons is being pursued 
using several approaches. Decay modes for which a CP- violating rate asymmetry 
would be particularly easy to interpret include those of neutral B mesons to CP 
eigenstates such as J/tp Kg and 7r"'"7r~. However, one is required to identify the 
flavor of the decaying meson: was it a B^{= hd) or B [— hd) at the time of 
production? 

One suggestion for neutral B meson identification involves the reaction 

e+ + e-^T(45)^S° + S^ , (1) 

in which the flavor of one meson is identifled by correlating it with the semileptonic 
decay products of the other. In this configuration, with C{BB) = — 1, any CP- 
violating asymmetry is odd in the difference of times of the two decays [1,2], so 
asymmetric storage rings may be needed to provide the necessary time resolution. 
It has also been proposed [2] to produce the BB pair at a slightly higher energy, 
in association with an extra photon from B* decay, in order to produce a state 
with C{BB) = 1 to overcome this difficulty. In both cases the leptonic "tag" of a 
neutral B gives rise to an associated pure B^ or !B°. 

Another method for determining the flavor of a neutral B relies on the semilep- 
tonic decay of a meson or baryon containing a h quark produced at high energies 
in association with the neutral B and with many other particles. Here the flavor 
of the decaying B is misidentified part of the time as a result of S° — mixing, 
but in contrast to the case of production at and just above the T(45') resonance, 
the initial state of the decaying neutral B is usually assumed to be an incoherent 
mixture of particle and antiparticle [3]. 

In order to avoid having to measure the decay of the associated particle, it has 
been proposed [4] to study the correlation of the decaying neutral B with charged 
pions produced nearby in phase space. This method has been used to identify 
neutral D mesons though the decays D*^ 7r~^D^ or D*~ — > n'D [5]. As in 
semileptonic tagging at high energy, the initial B^ and B^ states were assumed in 
Ref. [4] to be incoherent with one another. 

In the present paper we develop a way to identify a neutral B meson which 
embraces the above examples but allows one to test experimentally for arbitrary 
coherence properties of the initial state. The method relies upon the description 
of the initial tagged state in terms of a density matrix, in a two-component "qua- 
sispin" space spanned by B^ and B . The two-component description is standard 
for kaons (see, e.g., Refs. [6]), and has been applied as well to the case of 5's [7]. 

We have found a set of time-dependent measurements of decays of a tagged B 
which completely specifies the density matrix, permitting one to obtain phases in 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix through a study of CP-violating 
decay rate asymmetries. Each tagging procedure specifies a different density ma- 
trix. Previously it was necessary to anticipate the properties of such a matrix 
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with the help of theoretical assumptions. Our result avoids the need for any such 
procedure, and provides an answer to the question of coherence of the initial state. 

In brief, we advocate the measurement of B decays both to states of identifiable 
flavor (such as J/ipK*^, with K*^ K^tt^) and to eigenstates of CP (such as 
J/ihKs)- The obscrvables which specify the initial density matrix are: (1) a 
modulation amplitude of the time-dependent term; (2) a shift in phase of the 
time-dependent term, and (3) an overall intensity associated with the production 
of the CP eigenstatc. Once these quantities are measured, one not only determines 
the density matrix, but also obtains valuable information about weak phases. 

We now define the basis states and the density matrix. We choose 
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A density matrix p allows one to discuss incoherent and coherent states in a unified 
manner. The most general such 2x2 matrix has the form 

p^l[l + Q.^] , (3) 

where Q is a vector describing polarization in quasispin space, satisfying < 1, 
and ai {i — 1, 2, 3) are the Pauh matrices. We give two examples: 

(1) Any pure state corresponds to a linear combination of and B with 
arbitrary complex coefficients, whose sum of absolute squares equals unity. Such 
a state can be denoted by a density matrix with Q = |Q| = 1. 

(2) An arbitrary incoherent combination of B^ and 5^ with relative prob- 
abifities Pi and P2 = 1 — Pi corresponds to a diagonal density matrix with 
Q1 — Q2 — 0, Q3 — 2Pi — 1. This is the case we considered to hold in Ref. [4]. 

As a special case of either (1) or (2), one describes the density matrices for 
initial B^ and 5^ by diag(l,0) and diag(0,l), respectively. 

The probability for a transition from an initial state denoted by the density 
matrix to a final state denoted by p/ is then 

/(/) = TV (p^rV^T) , (4) 

where T is the amplitude which time-evolves the state from i to f. Here we shall 
take / to denote an arbitrary coherent superposition of S° and B at time t. We 
shall also be able to write a density matrix pf which takes account of the decay of 
this superposition. 

The timc-cvohition can be visualized by transforming from the B^, B basis to 
that of mass eigenstates. We shall neglect differences between lifetimes of these 
two eigenstates [1], and denote them by Bl ("fight") and Bh ("heavy"): 

\BL)^p\B') + q\^) , \BH)^p\B')-q\^) . (5) 
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Here p and q are approximately of unit magnitude [1]. We adopt a phase con- 
vention in which the S° — 5° mixing amphtude is real, so that p = q = 1/ V2- 
This differs from a more standard convention [1] by a phase which we shall take 
into account when calculating amplitudes for decays of b quarks. The transforma- 
tion between flavor eigenstates and mass eigenstates is then implemented by the 
unitary matrix 

^1 1 " 



V2 



^(^1 + ^3) . (6) 



1 -1 

The matrix describing the time evolution is diagonal in the mass eigenstate basis: 

^-iMnt ^ e-^*/Miag(e-^'"^*, e-^™«*) = e-r*/2g-imtgia3Amt/2 ^ 

where fh = {rriH + tol)/2, Am = niH — ttil. Thus the time evolution operator T 
in the S° basis is just T = U^e-'^^^^U = e-r*/2e-imtgiaiAmt/2_ 

The trace for the transition probability /(/) can be computed by applying the 
matrices U and W to the initial and final density matrices instead of to e~'-^°*. 
Defining p' = UpW, we find that the effect of U is to rotate Q into Q', where 

Q'i = Q3 , Q2 = -Q2 , Q's^Qi . (8) 

We already remarked that states which are pure B'^ or pure 5° correspond to 

Q2 — ±1, Qi = Q2 = 0. Their transformed density matrices are p'j = (1/2) (1 ± 
(Ti), respectively, since then Q[ = ±1, Q2 — Q'^ — 0. The transition probability 
can now be written in terms of traces as 

/(/) = Tr (p^e'^^'V/e-^^^*) . (9) 

The states and B^ at the time of decay t may be identified by their decays to 
states of identifiable flavor, e.g., 5° J/i/jK*'^, with K*^ K^-k' . We assume 
that a single weak subprocess contributes to the decay, which is an excellent 
approximation for these final states [1]. 

With the convention in which the mixing amplitudes p and q in the neu- 
tral B mass eigenstates are real, the weak decay amplitudes for B^ J/ipK*^ 
and W J/^K may be denoted Ae and Ae^^ , respectively, where /3 = 
Arg {—V*fycd/ytb^td): and Vij are elements of the Cabibbo-Kobayashi-Maskawa 
matrix specifying the charge-changing weak couplings of quarks. We find 

^|A|V^*[l±(g;cosAmi-g^sinAmi)] . (10) 

No reference to Q'^ appears. The sine and cosine terms may be combined into 
a cosine of a phase-shifted argument by defining 

Q[ = Q'^cosS , Q^ = Q1 sin 5 , (11) 
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so that 




^'A\^e^^'[l±Q'^cos{Amt + 5)] . (12) 



Whenever the initial state is any incoherent mixture of B'^ and B with relative 
probabilities Pi and P2 — 1 — Pi, respectively, one simply sets Q'j^ — 2Pi — 1 and 
S = in the above expression. 

We consider a charge- symmetric production process (such as pp or e+e^ col- 
lisions) in which an arbitrary coherent, partially coherent, or incoherent combi- 
nation of neutral B° and B is produced, with an additional particle of specific 
charge (such as a charged lepton or pion) bearing some specific kinematic relation 
to it. This relation could consist, for example, of a lepton with specific transverse 
momentum or range of transverse momenta, a definite-mass combination M{Bti) 
(or range of such combinations), a specific angle or range of angles between the B 
and tagging particle, or a specific rapidity difference (or range thereof) between 
the two particles. It is necessary to determine the relation of this process to the 
one in which the charge of the tagging particle is the opposite. 

Under the phase convention wc have chosen for B^ and B , if we take CP\Bl) — 
\Bl) and CP\Bh) — —\B}j), the charge conjugation operation has the phase 

C|5°) = -|#) ; C|5°) = -|5°) . (13) 

Under charge conjugation, the first and second rows and columns of the density 
matrix p are interchanged, so that Qi Qi, Q2 —Q2, Q3 —Qs- We then 
find that the effect of charge conjugation is a rotation by tt about the 3' direction: 

C: g; ^ -q;, ^ -g'^, g^ ^ g^ . (i4) 

Therefore, the decay rates /(/) for states tagged with antiparticles are then given 
in terms of those /(/) for states tagged with particles by 

Hf; Q[, Q2, Q3) = Hf; -Q'l, -Q2, Qs) ■ (is) 

For a final state identified as a S'^ by its decay to J/t/jK*^, the time-dependent 
asymmetry is 

As a consequence of the assumed charge symmetry of the production process, one 
has liJ/il^K*^) = liJ/i^K*^) and I{J/iliK*'^) = I{J /i^K*^). 

We can then measure both components Q'^ and Q^-, or, equivalently, Q'^ and 
5, using decays to fiavor eigenstates. The ALEPH Collaboration [8] has recently 
measured a time-dependent asymmetry of the above form, fitting it under the 
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assumption 6 = 0. It should be possible to set limits on 6 from this measurement. 
A way of measuring Q'^ will be described when we come to discuss CP-violating 
asymmetries. 

The above expressions are quite general and not restricted to any specific tag- 
ging method. The initial density matrix p can refer to any method of preparation, 
including specific configurations at or just above the T(45') resonance in e+e" 
collisions. 

A measurement of Q'^ for neutral nonstrange B mesons can be performed 
by utilizing their decays to the specific CP cigenstate J/ip Ks- In our phase 
convention, the amplitudes for and 5^ to decay into J/ipKs are A'e~^^/^/2 
and — A'e*^/ \/2. Here we have taken into account the intrinsic negative CP of 
the J/ipKs state, and neglected the small CP violation in the kaon system. The 
density matrix for the final state is then 



1/|2 



— e 



1/|2 



_g2i/3 



(17) 



For reference, we also note that 



PJ/i>KL 



2' ' 



-id 



2' ' 



,2i/3 



-2i/3 



(18) 



The expressions for the decay rates for states prepared with particle and an- 
tiparticle tags are then 




^|>l'|2e-^*{l± [Q3COs2/3 + Q'^sin2/3 sin(Ami + (^)]} , 



1 
2' 
1 
2' 



= ^|A?e-^*{l± [g3Cos2/?-g'^sin2/? sin(Ami + (5)]} . 



(19) 



As in the case of decays to the flavor eigenstates, the time-dependent term has 
a phase shift b and a modulation amplitude Q'^. It depends upon sin (Ami -|- 5) 
rather than cos(Ami -|- S). 

The decay asymmetry for the J/ i)Ks final state is then 



I{J/^Ks)-I{J/i)Ks) _ -Q'^ sin 2/? sin (Ami + 5) 



IiJ/iPKs)+IiJ/iPKs) 



1 - Q'3 cos 2/3 



(20) 



The component Q'^ (which appears even in the absence of CP violation) is a 
necessary ingredient in the discussion of possible coherence. It is this component 
that leads to correlations between Ks and Kl produced in (p decay, as discussed 
in Ref. [9]. In order to learn its value, we measure the rate for J/ipKs production 
(summing over particle and antiparticle tags, so that the time-dependent terms 
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cancel). We compare this with the corresponding sum of rates (which also has no 
time-dependence, and is independent of Q'^) for production of a flavor eigenstate 

Since we cannot observe a directly, we must resort to indirect means. If 
the rate of production of charged and neutral i?'s is equal, as is expected to be 
the case for high-energy e^e~ collisions [4], one can compare the rates for J /ipK^ 
and J/t/jKs production, making use of the fact that the decays of B mesons to 
J/ipK involve the quark subproccss h ccs, which conserves isospin [10]. 

A more general method makes use of the expectation that the K^/K^ ratio 
will be the same as the K*^/K*^ ratio. The flavor eigenstatcs involving K*^, K^, 
and K*~^ are all observable. Thus we can determine the ratio 

In either method we learn the relative normalization of rates for decays to 
flavor eigenstates and CP eigenstates. We thereby determine the magnitude of 
the term Q'^cos2f3, and then use the asymmetry in J/ipKs decay to measure 
sin 2/9. With the possibility of a discrete ambiguity (unlikely for known ranges of 
CKM parameters), we then obtain cos 2/3, thereby flnding Q'^ itself. 

The corresponding decay asymmetry for the tt'^tt^ flnal state is easily calcu- 
lated by our method. We neglect penguin effects [11], which can be dealt with by 
studying the 27r^ flnal state. With our phase convention for b quarks, the result 
can be obtained by the substitution f] — * —a in the corresponding result for the 
J/i/jKl flnal state, where a = Arg (-V^ftl^d/ICbKd)- We flnd 

, ^ _ /(tt+tt^) — /(tt+tt^) —Q'± sin2a sm{Amt + S) , . 

■>= /(tt+tt-) + /(tt+tt-) " l + g^cos2a ' ^ ' 

Since wc have already measured all components of Q' and the phase 5, this result 
can be used to extract a. 

To sum up, we have developed a density-matrix formalism for the discussion 
of decays of neutral B mesons which have been "tagged" in any arbitrary fashion. 
The density matrix is specifled in terms of a polarization vector Q in the basis 
labelled by B^ and B , or Q' in the basis of mass eigenstates. It is the com- 
ponents of this polarization vector which determine all observable CP-violating 
asymmetries, and which can be learned by measurements on the time-dependence 
of decays to J/ipK* and J/ipK. One learns as well the angle (3. Decays to fl- 
nal states such as 7r"^7r~ can then provide information on other phases of the 
Cabibbo-Kobayashi-Maskawa matrix. 
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